arXiv: 1505.02530v2 [cond-mat.mes-hall] 11 Oct 2015 


Current-induced spin polarization and spin-orbit torque in graphene 
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Using the Green function formalism we calculate a current-induced spin polarization of weakly magnetized 
graphene with Rashba spin-orbit interaction. In a general case, all components of the current-induced spin 
polarization are nonzero, contrary to the nonmagnetic limit, where the only nonvanishing component of spin 
polarization is that in the graphene plane and normal to electric field. When the induced spin polarization is 
exchange-coupled to the magnetization, it exerts a spin-orbit torque on the latter. Using the Green function 
method we have derived some analytical formulas for the spin polarization and also determined the correspond¬ 
ing spin-orbit torque components. The analytical results are compared with those obtained numerically. Vertex 
corrections due to scattering on randomly distributed impurities is also calculated and shown to enhance the spin 
polarization calculated in the bare bubble approximation. 

PACS numbers: 72.20.My, 72.80.Vp, 72.25.-b 


I. INTRODUCTION 

One of the main issues of the present-day spin electronics, 
that is of great importance for further development of high- 
density memory devices and magnetic random access mem¬ 
ories, is effective manipulation of magnetization by a spin- 
polarized current. The general idea of switching and control¬ 
ling orientation of a magnetic moment with electric current 
flowing through a system is based on coupling between the 
electron spin and magnetic moments. Two such interactions 
turned out to be especially useful - exchange interaction and 
spin-orbit coupling. 

In a magnetically nonuniform system, the spin-polarized 
current generates a torque that is a consequence of: (i) ex¬ 
change coupling between the conduction electrons and mag¬ 
netization, and (ii) conservation of angular momentum in the 
system. The torque appears then as a result of the spin an¬ 
gular momentum transfer from a spin-polarized current (or 
pure spin current) to magnetic moments. Therefore, this 
torque is called spin-transfer torqueSuch a torque leads, 
among others, to magnetic switching in spin vales and to do¬ 
main wall displacements, as observed recently in many exper¬ 
iments. Moreover, these phenomena give a possibility to con¬ 
struct low-power non-volatile memory cells (STT-MRAM, 
racetrack memory), integrated circuits employing a logic-in¬ 
memory architecture as well as logic schemes processing in¬ 
formation with spins 

Another possibility to control orientation of magnetic mo¬ 
ments is based on a spin torque that appears due to spin-orbit 
interaction in the system. The corresponding torque exerted 
on the magnetization is usually referred to as the spin-orbit 
torque, and appears also in a magnetically uniform system, 
like a single uniform layer. Physical mechanism of the spin- 
orbit torque is based on a nonequilibrium spin polarization 
of the system, which is induced by an external electric fleld 
(current) in the presence of spin-orbit interaction. Such a spin 
polarization was predicted long time ago in nonmagnetic sys¬ 
tems, where an electric current flowing through the system 
with spin-orbit interaction was shown to induce not only the 
transverse spin-current^^ (so-called spin Hall effect), but also 


a spin-polarization of conduction electrons In the case of 
two-dimensional electron gas with Rashba spin-orbit interac¬ 
tion, the induced spin polarization is in the plane of the elec¬ 
tron gas and normal to the electric fleld. Such a nonequilib¬ 
rium spin-polarization may be treated as an effective magnetic 
fleld, which may lead to reorientation of a magnetic moment, 
and also can modify or induce magnetic dynamics. The spin- 
orbit torque was analyzed in recent few years in many papers, 
mainly in metallic and semiconductor heterostructures 
While the current-induced spin polarization, known also as 
the inverse spin-galvanic effectr^ is well known and was in¬ 
vestigated theoretically as well as experimentally in the recent 
three decades, the role of geometric phase in this effect, and 
consequently in the spin-orbit torque, was invoked only very 
recentlyJ^-— 

In this paper we consider the current-induced spin polariza¬ 
tion and spin-orbit torque in graphene, which is assumed to 
be deposited on a substrate that ensures the presence of spin- 
orbit interaction of Rashba type^^. We also assume that the 
graphene is magnetized, which may be either due to the mag¬ 
netic proximity effect to a ferromagnetic substrate (or cover 
layer), or due to magnetic atoms (nanoparticles) on its sm- 
facei22rii Coexistence of the Rashba spin-orbit interaction and 
proximity-induced magnetism in graphene was predicted the¬ 
oretically and also observed experimentallyi2i-ii As the spin 
transfer torque in ferromagnetic graphene junctions was al¬ 
ready considered theoretically (see e.g. Yokoyama and Lin¬ 
der—), the problem of spin torques induced by spin-orbit in¬ 
teraction in graphene is rather unexplored. 

It has been shown that the current-induced spin polariza¬ 
tion in a defect-free nonmagnetic graphene with Rashba spin- 
orbit interaction is oriented in the graphene plane and is also 
normal to the current orientation. Moreover, sign of the spin- 
polarization depends on the chemical potential and also on 
the sign of the Rashba spin-orbit coupling parameteri^ When 
the Fermi level passes through the Dirac points, the spin po¬ 
larization becomes reversed. In this paper we show that the 
current-induced spin polarization in magnetic graphene has 
generally all three components. In the approximation linear 
with respect to the magnetization, one of these components is 
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equal to that in the case of a nonmagnetic graphene, i.e. it is 
proportional to the relaxation time. The leading terms in the 
other two components are independent of the relaxation time. 

The paper is organized as follows. In Sec. 2 we describe 
the model and present a general formula describing current- 
induced spin polarization. Analytical formulas as well as nu¬ 
merical results for the current-induced spin polarization are 
presented in Sec. 3. Vertex correction is calculated in section 
4, while the spin-orbit torque is described and discussed in 
Sec. 5. Summary and final conclusions are in Sec. 6 . 


II. MODEL AND METHOD 



FIG. 1. (Color online) Schematic of the system under consideration. 
Graphene is on a substrate which assures a nonzero magnetization 
and also a spin-orbit interaction of Rashba type. Orientation of the 
magnetic moment M is described by the angles 0 and An external 
electric field is oriented along the axis y. 


Transport properties of graphene close to the charge neu¬ 
trality point are determined mainly by electrons in the vicinity 
of Dirac points. The corresponding effective-mass Hamilto¬ 
nian, which describes the low-energy electronic states in 
graphene around the K point of the Brillouin zone, can be 
written as a sum of three terms ri 

= + + ( 1 ) 


The first term, i^o, describes the low energy electronic states 
of pristine graphene, and can be written as a matrix in the 
pseudospin (sublattice) space. 
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where v = Hvf, with vf denoting the electron velocity in 
graphene, which is constant. The second term in Eq. (1) de¬ 
scribes the Rashba spin-orbit interaction due to a substrate. 
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with A being the Rashba spin-orbit coupling parameter. The 
last term of the Hamiltonian (1) represents the infiuence of 
an effective exchange field M created by a nonzero magne¬ 
tization. Such a magnetization can appear in graphene, for 
instance, due to the proximity effect to a magnetic substrate. 
This term can be written in the form. 


its orientation will be described by two spherical angles, 0 and 
as indicated in Fig. 1. Moreover, the absolute magnitude of 
M is assumed to be constant, |M| = M = const. Hamil¬ 
tonian for the second non-equivalent Dirac point, K', can be 
obtained from Hr by reversing sign of the wavevector com¬ 
ponent kx and substitution ay —> —ay in Hr. 

In the lowest order with respect to the exchange field M, 
the casual Green function corresponding to the Hamiltonian 
( 1 ), = {[s ^ /j. + i5 sign ( 5 )] — has poles at 5 = 

En — — sign ( 5 ), where (n = 1 — 4) are eigenvalues of 
the Hamiltonian (1) without the term These eigenvalues 
have the following form: 

= TA - + ^2 (5) 

E3,4 = tA + (6) 

where £’ 1,2 correspond to the valance bands, while £ 3^4 de¬ 
scribe the conduction bands. Note, the bands corresponding to 
n = 2 and n = 3 touch each other at the Dirac point (k = 0), 
while a gap equal to 4A appears between the bands n = 1 and 

n = 4. 

In the presence of a dynamical external electric field applied 
along the axis y, the total Hamiltonian for electrons near the 
K point takes the form 

H = H% + H^, (7) 

where the second term. 
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where the exchange field M is measured in energy units. This 
field can be related to the magnetization M and the local ex¬ 
change interaction between the conduction electrons and mag¬ 
netization in the two-dimensional graphene, Jex(i* — r') = 
Jex<^(r — r'), via the formula M = (Jex/2^/iB)M, Here, g 
is the Lande factor (g = 2), fiR is the Bohr magneton, while 
positive and negative Jex (measured in the units of Jm^) cor¬ 
respond to antiferromagnetic and ferromagnetic coupling, re¬ 
spectively. In the above equations, cr is the vector of Pauli 
matrices, cr = {ax^ay^ az), while the matrix ao denotes the 
unit matrix in the spin space. Note that the so-called intrinsic 
spin-orbit interaction in graphene is very small and therefore 
it is neglected in our consideration. In a general case, the mag¬ 
netization vector M may be oriented arbitrarily in space, and 


is the perturbation due to interaction with the time-dependent 
electromagnetic field represented by the vector potential 
Ay{t) = Aye~'^^^. Here, e is the electron charge, Vy is the y- 
component of the electron velocity operator, v = OHr /5k, 
whereas uj is the frequency of the dynamical field (later we 
will take the limit of cc —> 0 ). 

When an electric current flows in the system due to the elec¬ 
tric field, electron spins become polarized as a result of the co¬ 
operation of the current and Rashba spin-orbit coupling. This 
nonequilibrium spin polarization of conduction electrons can 
be calculated (in the zero-temperature limit) using the follow¬ 
ing formula: 
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where t') is the zero-temperature causal Green function 
corresponding to the total Hamiltonian H (see Eq. 7), and Sa 
is the spin vertex function defined as 


S 
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h f (Ta 0 \ 

2 \ 0 (Ja ) ' 


( 10 ) 


Upon Fourier transformation with respect to the time variables 
and expansion in a series with respect to the vector potential 
Ay = —iEyjuj, the expression (0 for the induced nonequi¬ 
librium spin density takes the form 

S(y{uj) =J J —Sa G^{e E huj)vyG^{s).A^) 


In the dc limit, cc ^ 0, the above formula leads to the follow¬ 
ing expression for the spin polarization: 


Sa = —EyhTT 
Ztt 






( 12 ) 


where G^^^^ is the retarded (advanced) Green function cor¬ 
responding to the unperturbed Hamiltonian (1), taken at the 
Fermi level {e = 0). Upon taking into account Eqs ([TOb 
and (O, and also including the contribution from the second 
Dirac point, the expression for the induced spin polarization 
acquires the form 


eh^ 

-—EyTi 
27r ^ 




Gi'^VyGi^. (13) 


Based on this formula we calculate analytically as well as nu¬ 
merically the current-induced spin polarization, as described 
and discussed in the subsequent section. 


III. CURRENT-INDUCED SPIN POLARIZATION 


Up to the terms linear in the exchange field M, the func¬ 
tions Ta (a = z), see Eq. (15), can be written as follows: 


Tx = 16X7rvjii{k^v^ — 4A^/i^), (16a) 

Ty = 6A7rvXMy{k‘^v^ - 2A^)r sin 0, (16b) 

Tz = — 64:71 XMv^k^ juT cos6>sin^. (16c) 


Note, the dependence on the orientation of M is contained in 
the above expressions for Ty and T^, while Tx is independent 
of M. Equations (14) and (16) allow finding spin polarization 
in a general case, i.e. for an arbitrary relaxation time. How¬ 
ever, some analytical expressions for all components of the 
spin polarization can be obtained in the limit of low impurity 
concentration, i.e. for long relaxation times (r oo). 

Consider first the x-component of the spin polarization. 
C ombining E q. dlbab with Eq.([T^ and making the substitution 
y/k‘^v‘^ + A^ = 7 , one obtains 

Sx = SehEyXjii 

d-f-f (72-a2)2-;,4^4aV 

"" A t;(27r)2 [(m + A + 7 )^ + r2] + 7 )^ + P] 

^[(M + A-7)^ + r2][(;,-A-7)2 + r2 ]- 

From this formula follows that Sx is independent of M in the 
linear approximation with respect to the exchange field. For 
long relaxation times we get the same analytical formulas as 
those in the case of nonmagnetic graphene^^ i.e.. 




e 2A zb /i 


(18) 


for the Fermi level lying in the range — 2A < /i < 2A, and 


From the general formula ([T^ one finds the following ex¬ 
pression for the Qf-th component of the spin polarization: 
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where T„ is defined as 
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Here, is the nominator of the retarded (advanced) 

Green function, (j) stands for the angle between the axis x and 
the wavevector k, while F = h/2r, where r is the momentum 
relaxation time. The parameter F (or equivalently relaxation 
time r) will be treated here as a phenomenological parameter, 
which effectively includes contributions due to momentum re¬ 
laxation from various scattering processes (scattering on im¬ 
purities, other structural defects, phonons, or electron-electron 
scattering). Note that F depends in general on the chemical 
potential g and may be also different in the two Rashba sub¬ 
bands. However, when the Fermi level is in the two subbands, 
we assume for simplicity the same F for both of them. 


for |/r| > 2A. In both above equations (as well as below), 
the upper and lower signs correspond to /i > 0 and /i < 0, 
respectively. 

The spin polarization given by Eqs ([18]) and ([T9l) is propor¬ 
tional to r. However, one should bear in mind that these for¬ 
mulas were derived on the assumption of long r. Therefore, 
one may expect some deviations from this formula when r is 
finite and not too long. In Fig. [2a) we show variation of the 
Sx component of spin polarization with the chemical potential 
ji and relaxation time r, obtained by numerical integration of 
the formula (17). Figures IS) and|2|:, in turn, present cross- 
sections of the density plots shown in Fig. |2a) for constant 
values of r and (i, respectively. The results obtained from 
the analytical formulas are compared in Figs |2b) with those 
obtained by numerical integration of the formula (17). From 
this comparison follows that for r of the order of 10“^^s or 
smaller, there are some deviations from the results given by 
the analytical formulas, though these deviations are not large. 
For r of the order of 10“ ^^s or longer, numerical results match 
quite well those obtained from the analytical formulas. Since 
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FIG. 2. (Color online) Spin polarization components induced by a current: (a,b,c) show the x-component, (d,e,f) show the ^-component, 
whereas (g,h,i) show the ^-component. The top panel (a,d,g) shows the spin polarization components as a function of chemical potential ji 
and relaxation time r. The medium panel (b,e,h) shows the spin polarization components as a function of chemical potential /x for indicated 
values of r, while the bottom panel (c,f,i) shows the polarization components as a function of r for indicated values of /x. The right parts of 
(d,e,f) present the corresponding shaded regions in the left parts. The solid and dashed lines in (b) represent the results based on the analytical 
formulas and numerical integration, respectively. The curves for r ^ oo in (e,h) correspond to analytical solutions. The other parameters are: 
A = 2 meV, Ey = 1 V/cm, M = 0.1 meV, 0 = tt/S, and ^ = 7r/2. 


the Sx component is the same in magnetic and nonmagnetic 
limits (within the approximations used here), and in the non¬ 
magnetic limit it was considered and analyzed in Ref. (HI], we 
will not discuss this component in more detail. 

From Eqs (O and (I16bl) one finds the ^-component of the 
spin polarization in the following form: 

Sy = 32eEyhXiJ.M 

d-f-f (7^-3A2)rsin(9 

""A v{27rr [(Ai + A + 7)2 + - A + 7 )^ + n 


[(M + A-7)^ + r2][(;,-A-7)2 + r2] ■ ^ 

In the limit of slow relaxation, F -^0, the above formula 
leads to the following analytical results: 
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for \ii\ < 2A, and 


( 21 ) 
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for |/x| > 2A. 

Numerical results for the ^-component of the current- 
induced spin polarization, obtained by numerical integration 
of the formula (20) are shown in Fig.2(d) as a function of 
chemical potential /x and relaxation time r. Figures 2(e,f) 
present cross-sections of Fig.2(d). Figure 2(e) additionally 
shows the results obtained from analytical formulas, see the 
curves for r ^ oo. Right parts of Fig.2(d,e,f) present in more 
detail the corresponding shaded regions. Similarly as the x- 
component, Sy is antisymmetric with respect to reversal of 
the sign of Fermi energy, and its dependence on fi also re¬ 
veals some steps at /x = ±2A. These steps are associated 
with the edges of the bands Ei and E/^. Moreover, when the 
Fermi level is at the Dirac point (/x = 0), the analytical so¬ 
lution ([21]) for Sy becomes divergent. To understand origin 
of the divergency in the analytical solution for r ^ oo, one 
should note that the solution for the x-component is also infi¬ 
nite for r ^ oo, independently of /x. This clearly shows that 
the limit of r ^ oo is not physical as the dissipation processes 
are necessary in order to stabilize a finite current-induced de¬ 
viation of the system from equilibrium, and thus also a finite 
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current density and spin polarization. Therefore, in Fig.2(e) 
we compare the numerical results based on the corresponding 
analytical formulas with those obtained by numerical integra¬ 
tion. This comparison clearly shows that the results obtained 
from the analytical formulas are roughly in agreement with 
those obtained from numerical integration, except the vicinity 
of /i = 0, where the analytical solution diverges for /i ^ 0, 
while the numerical results based on Eq.(l20t are then finite. 
Moreover, some discrepancy also occurs around /i = ±2A, 
but now the difference is finite and rather small. Thus, one 
should bear in mind that the analytical results ([21]) and (l22l) 
for the ^-component have limited applicability range, and are 
not applicable for /i in the vicinity of the Dirac points. 

The Sz component can by found from Eqs (O and (I16cl) 
and acquires the form 


The other two components behave more regularly in this re¬ 
gion. All three components exhibit some cusps (or dips) when 
fi is in the vicinity of/i = ±2A,i.e., when the Fermi level ap¬ 
proaches the top edge of the band Ei or bottom edge of the 
band E^. The spin polarization also remarkably depends on 
the Rashba parameter A. This dependence reveals peculiari¬ 
ties of the corresponding electronic structure, and remarkably 
depends on the Rashba parameter. In numerical calculations 
we assumed the Rashba spin-orbit coupling parameter A = 2 
meV. Generally, this parameter depends on the substrate (or 
cover layer), and in real systems varies from a few to a few 
tens of meV, see eg. Refs Isslfioll . 


IV. VERTEX CORRECTION 


Sz = —32eEyhXiJ.M cos 0 sin ^ 

^ r Jen _ (7^ - A^)r _ 

J\ i'(27r)2 [(^ + A + 7)2 + r2] [(// - A + 7)2 + r2] 




Similar calculations as those done for the ^-component lead 
to the following analytical expressions in the limit of long re¬ 
laxation time: 


= (24, 

for ji < 2A, and 




ehM . . , 

T";-^ cosa'sm^ 

Att a 


11 ^ - 2A^ 
— A^) 


Eij 


(25) 


for ji > 2A. 

In Fig.2(g) we present the 2 )-component of the current- 
induced spin polarization, calculated as a function of the 
chemical potential and relaxation time by numerical integra¬ 
tion of the formula (25). In turn, Figs.2(h,i) show cross- 
section of Fig.2(g). In Fig.2(h) we additionally compare the 
numerical results with those obtained from analytical solu¬ 
tion. Now, the analytical solution is not divergent, see the 
curve for r —> oo. When the relaxation time is sufficiently 
small, the numerical results obtained from Eq. (1231) deviate 
from the results obtained on the basis of the analytical formu¬ 
las. These deviations are rather small forr > 10“^^s, except 
the region near the zero chemical potential. However, the dif¬ 
ference between the analytical and numerical results around 
/i = 0 is now much less pronounced than it was in the case of 
the ^-component (compare Fig.2(e) and Fig.2(h). In turn, for 
r < 10“^^sthe deviations become remarkable in the whole 
range of the chemical potentials shown in Fig.2(h). 

All the components of the spin polarization (Sx, Sy and 
Sz) vanish at /i = 0 and are antisymmetric with respect to 
the sign reversal of the chemical potential. Numerical results 
presented above show that the spin polarization strongly de¬ 
pends on the Fermi level position. In the close vicinity of 
the Dirac points, the ^-component of the spin polarization has 
pronounced peaks (positive above and negative below /i = 0). 


In the preceding section we have calculated spin polariza¬ 
tion induced by electric field assuming effective relaxation 
time r (or equivalently relaxation rate F). Both, r and chem¬ 
ical potential were treated there as independent parameters. 
When considering a specific relaxation mechanism, these pa¬ 
rameters usually are not independent. Since the dominant 
scattering processes are on impurities, we consider now this 
problem in more details. Assume the scattering potential 
created by randomly distributed weak short-range scatterers, 
which may be written as V (r)socro with Gaussian correlations 
(IA(r)I/(r'))) = niV‘^6{T — r') (where sq and ao and de¬ 
note unit matrix in the pseudo-spin and spin subspace respec¬ 
tively). 

Detailed calculation of the self energy due to scattering on 
the point-like impurities gives ri ^4 = — A) and 

r 2,3 = ^ 2^2 (ImI + A), where rii is the impurity concentration 
while V is the impurity scattering potential. When |/i| ^ A, 
then indeed ri ,4 r 2,3 = F. Otherwise, we take F as the 
average of ri,4 and r 2 , 3 , i-e. F = 

When calculating the impurity averaged conductivity, it is 
well known that non-crossing diagrams give an important con¬ 
tribution and renormalize the results obtained in the bare bub¬ 
ble approximation. Such a vertex renormalization is known 
to have a significant influence on the spin current induced 
via the spin Hall effect. In the case of two-dimensional 
electron gas with Rashba spin-orbit interaction it totally can¬ 
cels the spin Hall conductivity obtained in the bare bubble 
approximation.—'^ However, this is not a general property 
and in other systems the vertex corrections can only reduce 
partly the spin Hall effect.—"— 

The problem of disorder in graphene was discussed in many 
papers.—"— However, there is still a lack of information on the 
influence of disorder and impurities on spin-orbit driven phe¬ 
nomena in graphene. This problem was raised by Sinitsyn et 
al^ and Gusynin et alM in the context of spin Hall and spin 
Nernst effect in the presence of intrinsic spin-orbit interaction 
in graphene and in the case of spin-independent random po¬ 
tential. In this case problem becomes simpler because one 
can reduce the model to 2 x 2 space. Such a simplification, 
however, is not possible in the presence of Rashba spin-orbit 
interaction. 
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In the weak scattering limit, the localization corrections are 
vanishingly small and therefore only noncrossing ladder dia¬ 
grams are important. The summation over the ladder diagrams 
can be represented by the vertex corrections to the current- 
induced spin polarization. The renormalized spin vertex func¬ 
tion is then given by the following equation:— 


In turn, for a = 2 ; we find 

az = bz = = 0, 

^ 1 - riiV^Iz ^ 

where 


(35) 

(36) 


Sa — Sa + TliV' 


/ 




G A q 
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dkk _ Xz(A^,r) _ 
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where Sa is defined by Eq. (10). For the point-like scattering 
potential one can postulate the vertex function Sa in the form 


Sa — 


+ Cq 


h 

' 2 
h 



0 

0 


+ be 




(27) 


for a = x^y^z, where ba, Ca, and da are certain param¬ 
eters to be determined. To find these parameters we multiply 
Eq. (26) by the matrix as specified below and take the trace. 



for i = 0, X, y, z. Taking into account Eq. (27), one finds then 
a set of equations for the coefficients Ra^ba^Ca^da- 

We recall that in this paper the exchange field due to prox¬ 
imity effect is assumed to be small, so the current-induced 
spin polarization is limited to the terms linear in the exchange 
field. Consequently, the vertex correction is also calculated in 
the lowest order appropriate to have spin polarization linear in 
M. 

For a = X we find that: 


where 


bx — — dx — 0, 

1 


1 - niV‘^Ix ’ 


(29) 

(30) 


T 

-^x 


/ 


dkk _ XxifJ-X) _ 

27r nLi(Ai -En + -Er,- iT) 


(31) 


and 


Xx{^^, r) = + kA^{3T^ - ^ 2 ) + (r 2 + ^2)3 

+4(r^ - /)a2 + k^v^{r^ + //2)(3r2 + 4X‘^) 

Si (fc2i;2 _ ^2^ 4^2^2^ 

+ {3k\^ + 3 / + 2k^v^ { 1 ? + 2X^)) r^.(32) 

In the above equation, only terms up to the second order in 
r have been retained, while terms of higher order have been 
omitted. 

For a = y we find the following coefficients: 

Ciy = Cy = dy = 0, (33) 

by=az;= r). (34) 


and 


r) = (3r2 - ^ 2 ) (r 2 + ^ 2 ) 

+ (r^ + + k\'^ (sr^ - - 2X^) 

+ (r^ + (2 (r^ - 3/Li2) + 8A^) 

« ((fc2t,2 _ ^2)2 _ 4^2;^2)) ^^2^2 ^ ^2 _ 

+ {3k‘^v'^ + 2k^v^i? + 3/ - + 8A^) T^. (38) 

Finally the renormalized spin-vertex functions are: 



(39) 


(40) 


(41) 


This means that the results obtained in the bare bubble ap¬ 
proximation should be multiplied only by a numerical factor 
to take into account the vertex corrections due to disorder. 
More specifically, the results for Sx and Sy should be mul¬ 
tiplied by the factor y while those for Sz should be multiplied 
by (. The situation is significantly different from that found in 
the case of spin Hall effect. This is because transport phenom¬ 
ena and spin polarization are affected by scattering on impu¬ 
rities in remarkably different ways. 

In Fig.3(a) we show the renormalization parameter 77 as 
a function of chemical potential and relaxation time. Now 
the relaxation time is connected with the chemical potential 
through the relation ^ |/i|. A single point in the r, y 

space corresponds to a well defined value of niV‘^. How¬ 
ever, possible values of niV‘^ have been limited in Fig.3 to 
< (rziy^)max, where (nil/^)niax is a certain maximum 
value which is physically reasonable. The central white re¬ 
gion is bounded by the condition h/r = {niV‘^)ma.^\y\/x‘^ 
and is excluded for the considered parameters. In Fig.3c, in 
turn, we show the parameter 77 as a function of the relaxation 
time and the ratio niV‘^/{niV‘^)ena.x^- As one might expect, 
this figure shows that the normalization parameter 77 becomes 
reduced with decreasing niV‘^. Figures 3(b) and 3(d) present 
cross-sections of Fig.3(a) and 3(c), respectively. The above 
described results for 77 show that the Sx and Sy components 
are remarkably renormalized by the vertex correction and are 
enhanced by a factor of the order of 2 (between 1 and 3). 









7 


(a) 


12 

10 

8 




-20 


(b) 


-10 0 10 
/i [meV] 


20 



\li\ [meV] 




nyV(ny^)max 



1.6 
1.5 
1.4 
^ 1.3 
1.2 
1.1 
1.0 

2 4 6 8 10 12 14 161820 

ImI [iieV] 



c 




N 


0.01 0.25 0.5 0.75 


(h) 


nyV(ny2)max 


1.0 



nyV(ny^)max 


FIG. 3. (Color online) The parameter 77 as a function of the chemical potential and relaxation time (a) and as a function of relaxation time and 
riiV^/ (c). Figures (b) and (d) show 77 as a function of chemical potential /x (b) and riiV^/ (d) for indicated values 
the relaxation time. Figures (e)-(h) show the same variations as figures (a)-(d), but for the parameter The white regions in (a) and (e) are 
excluded for the assumed value of {niV^)ma^ = 0.4 x 10“^ (eV • nm)^. The other parameters are as in Fig.2. 


This enhancement of the spin polarization is comparable to 
that found in the case of two-dimensional electron gas with 
Rashbe interactioni^ The parameter in turn, is shown in 
Fig.3e-f. It is of the same order of magnitude as the parameter 
77 and depends on the chemical potential and relaxation time 
in a similar way, so we will not discuss it in more detail. 


V. SPIN-ORBIT TORQUE 

The current-induced spin polarization is exchange-coupled 
to the local magnetization M and thus exerts a torque on M. 
According to Eq.(|4|), energy of this interaction per unit area 
can be written as E’ex = — ( 2 /fi)M • S, where S is the induced 
spin polarization. Taking into account the relation between M 
and M, one finds the spin-orbit torque per unit area, r, exerted 
on the magnetization (more precisely on the corresponding 
equilibrium spin polarization of the system) in the form 

T = |m X S = X S. (42) 

h gi^B^ 

Let us consider in more detail some specific situations as 
concerns relative orientation of the magnetization and elec¬ 
tric field (current). Let us start with the situation when the 
magnetization M is in the plane of the system and perpendic¬ 
ular to the current. This corresponds to 0 = 0 and ^ = 0 
(Mx = M ^ 0 and My = Mz = 0). Lrom the above general 
equation follows that the spin-orbit torque can be then written 


in a general form as 

T = A{-jMxSz + kMxSy), (43) 

where i, j, and k are unit vectors along the axes x, y and 
z, respectively, and we introduced the following abbreviation: 
^ = Je^/ggsk. Taking into account Lqs (1211) . ([22b . ([24b and 
([25b . one finds immediately that the spin-orbit torque in this 
geometry disappears because both Sy and Sz component of 
the spin polarization vanish. 

Consider now the situation corresponding to 6 > = 0 and 
^ = 7 r /2 (My = M 0 and Mx = Mz = 0), i.e. the 
case when the magnetization is parallel to the electric current. 
Lrom Lq.([42b follows that the spin-orbit torque has the fol¬ 
lowing general form: 

T = A(iMySz — kMySx)- (44) 

The Sz component is now nonzero, and thus both, Sz and Sx 
contribute to the torque in this geometry. 

When = 7 r/ 2 , namely Mz = M 7 ^ 0 and Mx = My = 0, 
the magnetization is perpendicular to the graphene plane. The 
spin-orbit torque takes then the general form, 

r = A(-iMzSy^jMzSx). (45) 

Similarly as in the preceding situation, both Sy and Sx are 
nonzero and determine the torque. 

In the last two cases the spin-orbit torque contains two com¬ 
ponents: linear term with respect to Jgx (proportional to Sx) 
and quadratic term in Jex (proportional to Sz and Sy). The 
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spin orbit torque contains one component proportional to the 
relaxation time and another component whose the dominant 
part is independent on the relaxation time. 

In a general case of arbitrary orientation of the magnetic 
moment, magnitude and character of the spin-orbit torque 
varies with the orientation of the magnetic moment. This is 
because two components of the current-induced spin polar¬ 
ization depend on the magnetization, while the third one is 
independent of M. As a result the spin torque may have field¬ 
like and (anti)damping terms. 

VI. SUMMARY 

We have calculated current-induced spin polarization in 
graphene deposited on a ferromagnetic substrate, that ensures 
not only Rashba spin-orbit interaction but also a ferromag¬ 
netic moment in the graphene layer. To describe electronic 
spectrum of graphene we have used Kane Hamiltonian that 
describes low-energy states around the Dirac points. Using 
the zero-temperature Green functions formalism and linear 
response theory, we have derived analytical formulas for the 
spin polarization, up to the terms linear in M. Numerical re¬ 
sults based on the analytical formulas have been compared 
with those obtained by numerical integration procedure. From 
this comparison we have formulated applicability conditions 


of the analytical results. Significant deviations of the analyt¬ 
ical results from those based on numerical integration have 
been found for relaxation times smaller than 10“^^s. 

The nonequilibrium (current-induced) spin polarization ex¬ 
erts a torque on the magnetization via the exchange interac¬ 
tion. This torque contains a term which is proportional to the 
x-component of the induced spin polarization and therefore 
is proportional to the momentum relaxation time. The torque 
also includes a component whose main part is independent of 
the relaxation time. 

The spin-orbit torque due to the interplay of external 
electric field and Rashba coupling at the interface between 
graphene and a magnetic layer can be used for instance to trig¬ 
ger magnetic dynamics and/or magnetic switching. Indeed, 
such a switching was observed experimentally in a recent pa¬ 
per by Wang et al^ However, instead of graphene they used 
M 0 S 2 - another two-dimensional honeycomb crystal. 
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